Branching ratio measurement of Ks 77 
decay using a pure Ks beam in the KLOE 

detector 



The KLOE Collaboration 

F. Ambrosino^, A. Antonelli^, M. Antonelli^, F. Archilli^, 
C. Bacci^, P. Beltrame^, G. Bencivenni S. Bertolucci^, 
C. Binis, C. Bloise% S. Bocchetta^, V. Boccis, F. Bossi% 
P. Branchini^, R. Caloi^, P. Campana^, G. Capon ^, 
T. Capussela^ F. Ceradinii^, S. Clli^ G. Chiefari^ 
P. Ciambrone^, E. De Lucia ^, A. De Santis®, P. De Simone^, 
G. De Zorzis, A. Denig^, A. Di DomenicoS, C. Di Donato^, 
S. Di Falco^ B. Di Micco^, A. Doria^, M. Dreucci% G. Felici% 

A. Ferrari^, M. L. Ferrer^, G. Finocchiaro ^, S. Fiore^, 
C. Forti^, P. Franzinis, C. Gatti^, P. Gauzzi^, S. Giovannella ^, 
E. Gorini*^, E. Graziani^, M. Incagli*^, W. Kluge^, V. Kulikov'', 

F. Lacava^, G. Lanfranchi^, J. Lee-Franzini D. Leone ^, 
M. Martini ^ P. Massarotti^, W. Mei^, S. Meola^, S. Miscetti^ 
M. Moulson^, S. Miiller^, F. Murtas^, M. Napolitano^ 
F. Nguyen'^, M. Palutan^, E. Pasqualucci A. Passeri'^, 
V. Patera ^'^ F. Perfetto^, M. Primavera P. Santangelo^, 
G. Saracino^, B. Sciascia^, A. Sciubba^'^, F. Scuri^, L Sfiligoi^, 
T. Spadaro% M. Testa §, L. Tortora^, P. Valente^, 
B. Valeriani^, G. Venanzoni^, R.Versaci^, G. Xu^'j 

^Laboratori Nazionali di Frascati delVINFN, Frascati, Italy. 

^Institut fiir Experimentelle Kernphysik, Universitdt Karlsruhe, Germany. 

^ Dipartimento di Fisica dell'Universitd e Sezione INFN, Lecce, Italy. 

Dipartimento di Scienze Fisiche dell'Universitd "Federico IF' e Sezione INFN, 

Napoli, Italy 

^ Dipartimento di Fisica delVUniversitd e Sezione INFN, Pisa, Italy. 

^ Dipartimento di Energetica dell'Universitd "La Sapienza" , Roma, Italy. 

^ Dipartimento di Fisica dell'Universitd "La Sapienza" e Sezione INFN, Roma, 

Italy. 



^ Dipartimento di Fisica deU'Universita "Roma Tre" e Sezione INFN, Roma, Raly. 

^Physics Department, State University of New York at Stony Brook, USA. 

^Permanent address: Institute of High Energy Physics of Academica Sinica, 

Beijing, China. 

^Permanent address: Institute for Theoretical and Experimental Physics, Moscow, 

Russia. 



Abstract 

We have analyzed 1.62 fb"^ of e+e" collisions at a center of mass energy ~ 
collected by the KLOE experiment at DA^NE. This sample corresponds to a pro- 
duction of ~ 1.7 billion of Kg Kl pairs which allowed us to search for the rare 
~^ 77 decay. Kg are tagged by the Ki^ interaction in the calorimeter and the 
signal is searched for by requiring two additional prompt photons. Strong kine- 
matic requirements reduce the initial 0.5x10^ events to 2300 candidates from which 
we extract a signal of 600 it 35 events. By normalizing to the Ks 2tt^ de- 
cays counted in the same sample, the measured value of BK{Ks — > 77) is (2.27 
±0.13{stat.)toM{syst.)) x 10"^, in agreement with ©(P"') Chiral Perturbation The- 
ory predictions. 

Key words: e'^e" collisions, DA$NE, KLOE, rare Kg decays, xPT 



1 Introduction 

A precise measurement of the Ks — > 77 decay rate is an important test of 
Chiral Perturbation Theory (xPT) predictions. The decay amplitude of Ks 
'-f'-f has been evaluated at the leading order of xPT [1], 0{P^). providing the 
estimate of the corresponding branching ratio, BR[KS 77) = 2.1 x 10^^, 
with a few percent precision. This result is in agreement with the experimental 
measurement of NA31, that obtained BR{Ks 77) = (2.4 ± 0.9) x 10"^ 
[2]. The last precise determination of this BR{Ks '~f'~f) of 2.71 x 10^^, 
with a total uncertainty below 3%, comes from NA48 [3]. The last mentioned 
result differs from xPT 0{P'^) prediction of about 30%, indicating possible 
contributions from higher order corrections. 

In this paper, we show our search based on a data sample of 1.62 fb~^ of 
e+e- collisions collected with the KLOE detector [4] - [7] at DA$NE [8], the 
Prascati 0-factory. DA$NE is an e~^e~ collider which operates at a center 
of mass energy, W, of ~ 1020 MeV, the mass of the 0-meson. Equal-energy 
positron and electron beams collide at an angle of (tt - 25 mrad) producing 
0-mesons nearly at rest. 0-mesons decay 34% of the time into nearly coUinear 
K^K pairs. Since J^^(0) — 1 , these pairs are in an antisymmetric state 
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so the final state is always KsK^. All these imply that detection of a Kl 
guarantees the presence of a Ks of given momentum and direction. KLOE 
takes advantage of this to identify iir^-mesons independent of the decay mode. 
We refer to it as Ks tagging. 

The data sample analyzed corresponds to a production of ~ 1.7 billions of 
KsKl pairs. In the analysis an equivalent statistics of simulated events for the 
background was produced, as well as a sample with simulation of the signal 
with a factor ~ 15 larger. Using these samples allow us to reach a statistical 
error of 5.6% on the signal. While this accuracy is statistically inferior to the 
most precise NA48 result, this new measurement, having completely different 
background composition and origin of systematics, as well as from a pure Ks 
beam, can help to clarify whether 0{P^) contribution are present. 



2 The KLOE detector 



The KLOE detector consists of a large cylindrical drift chamber, DC [4], of 
4 m diameter and 3.3 m length with an helium-based gas mixture, surrounded 
by a lead-scintillating fiber electromagnetic calorimeter, EMC [5]. A supercon- 
ducting coil around the EMC provides a 0.52 T field. Permanent quadrupoles 
for beam focusing are inside the apparatus and are surrounded by two com- 
pact tile calorimeters with veto purposes, QCAL [6] . In this analysis, only the 
calorimeter system is used. 

The calorimeter is divided into a barrel and two endcaps covering 98% of the 
solid angle. The modules are read out at both ends by photomultipliers with 
a ~ 4.4x4.4 cm^ readout granularity, for a total of 2440 cells. Both amplitude 
and time signals are time signals collected at the two ends. The amplitude gives 
a measure of the energy deposited in the modules and the time signal yields 
both the arrival time of particles and the position in three dimensions of the 
energy deposits. Cells close in time and space are grouped into a "calorimeter 
cluster". The "cluster energy" E is the sum of the cell energies. The cluster 
time T and position R are energy-weighted averages. Energy and time res- 
olutions are aE/E = 5.7%/-^E (GeV) and at = 57 ps/-^E (GeV) © 50 ps 
respectively. 

The QCAL detector comprises two tile calorimeters of ~ 5X° placed thick- 
ness close to the IP, surrounding the focusing quadrupoles. Each calorimeter 
consists of a sampling structure of lead and scintillator tiles arranged in 16 
azimuthal sectors. The readout is done via wavelength shifter (WLS) fibers 
coupled to mesh photomultipliers. The special arrangement of WLS fibers al- 
lows also the measurement of the longitudinal coordinate by time differences. 
The tiles are assembled to maximize efficiency for photons coming from the 
Kl decays, yet retaining a high efficiency also for photons coming from the 
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IP. The QCAL solid angle coverage is 0.94 < \cose\ < 0.99. 

The standard KLOE trigger [7] uses calorimeter and chamber information. For 
this analysis only the calorimeter signals are relevant. Two energy deposits 
with E > 50 MeV for the barrel and E > 150 MeV for the endcaps are 
required. Identification and rejection of cosmic-ray events are also done at the 
trigger level. 



3 Search of Ks — > 77 with a pure Kg beam 

3.1 Ks tagging 

At the center of mass energy of M^f,, the mean decay lengths of the Ks and Ki 
are ~ 0.6 cm and Xl ~ 340 cm respectively. About 50 % of K^s reach the 
calorimeter before decaying. -fC^'s are tagged with high efficiency (~ 30%) by 
identifying a Kl interaction, which we call " Xi-crash" . This " i^^^-crash" has a 
very distinctive signature in the calorimeter, given by a late {^k = 0.2) high- 
energy cluster un-associated to any track. The "i^^-crash" provides a clean 
Ks tag. In this analysis the fake-tag contribution is essentially negligible. The 
average value of the center of mass energy W, is obtained with a precision of 30 
keV for each 100 nb"^ running period, by reconstructing large angle Bhabha 
scattering events. The value of W and the "/Ti-crash" cluster position allow us 
to establish, for each event, the trajectory of the Ks with an angular resolution 
of 1° and a momentum resolution of ~ 2 MeV. 

We use for normalization the measurement of the dominant neutral decay 
mode of the Ks 27r°, always tagged by "A'^-crash". 

In the analyzed sample, we have ~ 480 xlO® Ks tagged events. Using the 
most precise value of BRKS —>■ 77, we expect ~ 1300 Ks —>■ 77 events to 
be produced and tagged. The advantage of the present measurement is that, 
by tagging, we can completely neglect the Kl — > 77 background, which is the 
major contamination in NA48 analysis. 

3.2 Simulation of background and signal 

The main expected background in this search are Ks 27r° events with two 
lost photons. Such losses can be due to photons 1) either out of acceptance, not 
reconstructed by the calorimeter or 2) merged together. For the simulation of 
the background Monte Carlo, MC, we use a production of KsKl decays 
corresponding to an equivalent statistics of ~ 1.1 fb~^. For the MC signal we 
used instead, a production equivalent to ~ 18 fb~^. 
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Fig. 1. Angular distribution of the generated photons for Ks — > 271*^ events af- 
ter "JCi-crash" tagging: photons reconstructed by the EMC calorimeter (top), lost 
photons (bottom). 

The photon properties in the simulation (resolutions and detection-efficiency) 
have been tuned with data using a large sample of tagged photons in — >■ 
7r"'"7r~7r° events [9] selected using only drift chamber information. The presence 
of additional clusters in the events due to accidental overlap with machine 
background, dafbkg, clusters has been taken into account by inserting these 
events in the MC at hit level in the simulation. The dafbkg insertion includes 
its rate dependence along the running period, and all other basic running 
conditions, such as beam parameters (W, P^, X^jj) which vary run by run. The 
interaction of the Kl in the calorimeter is also properly simulated. 

3.3 Event preselection 

After tagging, we select events by counting the number of prompt photons, 
N^, i.e. neutral clusters in the EMC having a time of flight consistent with 
a particle with (3 = 1 coming from the interaction point. For signal counting 
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wc require = 2. while for normalization purposes we require iV^ = 4. 
A tight constraint on f3 is imposed to reduce the effect of event losses due to 
accidental clusters from dafbkg. Moreover, to improve the rejection of the main 
background, Ks — > 27r° with two lost photons, we accept clusters with energy 
above 7 MeV and produced in a large angular acceptance, |cos(^)| < 0.93. 
After this acceptance selection, the distribution of the lost photons for the 
background is peaked in the forward direction, as shown by the simulation in 
Fig. 1. 

To improve the background rejection we require a veto from the QCAL calorime- 
ter. This veto consists of rejecting events having at least one hit in QCAL with 
energy above the pedestal and in time with collisions. In Fig. 2, the data dis- 
tribution of the difference, ATq, between the reconstructed time of the QCAL 
hits, Tq, and their expected time of flight, TOFq, is shown for all tagged 
events with A^^ = 2. The flat distribution of hits due to machine background 
events shows clearly separated peaks bunched with the RF period. The sharp 
in-time peak observed is instead due to the lost photons from the process 
Ks 27r° impinging on QCAL. We veto all events in a time window, TW, 
deflned as |ATq| < 5 ns. The QCAL veto successfully rejects ~ 70 % of the 
background while retaining a very high efficiency (~ 99.96 %) on the signal. 

Since we have not simulated the dafbkg events in QCAL , when applying the 
veto on data we should correct for the losses due to the accidental coinci- 
dence of these background hits in the time window used. A data-calibrated 
correction, Cq — 1 — Pq^-, has been developed by determining -Pq^ i.e. the 
probability to find a spurious hit in TW. To estimate it we use two out-of- 
time windows one before, ear/y, and one after, /ate, the collision time. The 
average value of the probability in these two control windows provides a first 
evaluation of the correction. Furthermore, to assign a systematic error to this 
determination, we have measured these losses also in a control sample with 
a "i^'i-crash" and a well reconstructed Ks n~^n~ decay. This last sample 
does not have any photons impinging on QCAL, thus allowing us to calculate 
directly the losses in TW. The ATq distribution for these events is almost 
flat as shown by the overlapped distribution (points) in Fig. 2. We flnally 
determine this probability to be: Pj^ = (3.51 ± 0.04,^^^ ± 0.26syst.)%- 

At the end of the acceptance and QCAL veto selection, we count 157 xlO^ 
events in data. 



3.4 Kinematic fitting and event counting 

To improve the signal over background ratio (S/B), we apply a kinematic flt 
procedure which imposes the conservation of the 4-momentum at the origin, 
the Ks mass and /? = 1 for each photon (Ndof=7). The flt uses the knowledge 
of the Kl 4-momentum provided by the "Xi-crash" position, W and P^. In 
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Fig. 2. Inclusive distribution of the difference between the measured arrival time and 
the expected time of flight of prompt photons in QCAL. All hits for events tagged 
by a Ki-cvash with = 2 (black solid line) or with a reconstructed Ks vr+vr^ 
decay (points). 

Figs. 3. left (.right) the distribution of the of this procedure, Xfit^ is shown 
for data and MC after acceptance selection without (with) the application of 
the QCAL veto. A large amount of the background is distributed at high Xfit 
values while the signal is contained in the low Xfit region. In the following, 
we retain events by cutting at Xfit < 20. The Monte Carlo estimates that the 
signal (background) efficiency of this cut is 63.3 % (0.5%). The S/B greatly 
improves, from 1/500 to 1/3. 

Two other variables with a high discriminating power against background are 
the invariant mass of the two photons, M^^, and the opening angle between 
the two photons in the Ks center of mass system, 6*^. Since the kinematic 
fit imposes the Ks Kl direction and the Ks mass, we use the reconstructed 
variables before fit constraining. In Fig. 4.top(bottom) the 2D-plot of M^^ as 
a function of cos{6*^) is shown for background (signal) events. In Fig. 5, the 
distribution of the same scatter-plot for data is also shown. 
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Fig. 3. Distributions of Xfit events after tagging and = 2 requirement: 

before (left) and after (right) the apphcation of the QCAL Veto. 

Before fitting the data with MC shapes, we have tested our simulation abihty 
of reproducing the signal by comparing with data a control sample of Kl 77 
decaying near the beam pipe and tagged by Ks vr+vr" events. Around 200 
pb~^ of data and 450 pb~^ of Monte Carlo have been used. A kinematic fit- 
procedure, similar to the one of the Ks —>■ 77 sample, has been applied. The 
background is reduced to a negligible quantity by retaining only the events 
with XfitO^^) < 20. a gaussian fit to the M^^ distribution of this control 
sample provides a central value of (496.2 ± 0.8) MeV in data and of (488.7 ± 
0.5) MeV in Monte Carlo. This corresponds to an average energy-scale shift 
of ~ 1% in the simulation. This motivated an in-depth data-MC comparison 
of energy response and resolution as a function of the incident photon energy. 
This study has been carried out by looking at the energy pulls of the kinematic 
fit with a sample of ~ 80 pb~^ of Ks — >■ 27r° spread out over the entire data 
taking period. An ad-hoc correction has been applied to better calibrate the 
simulation versus the data. After applying this correction, the comparison 
between data and MC for the Kl — >■ 77 control sample is improved as shown 
by the fit to the M^^ distributions shown in Fig. 6. This control sample has 
been also used to make a data-MC comparison for the XfitO^^) distribution. 
A good agreement is also observed for this variable. 

For the i^'s'-tagged events, a binned maximum-likelihood fit to the 2D cos{9*^), 
M^^ distribution on data is performed relying on the MC signal and back- 
ground shapes. The likelihood function properly takes into account data and 
MC statistics of the distributions. The resulting x^/Ndof of the fit is 1.2. 
The quality of the fit is shown in Fig. 7 by comparing with data the simu- 
lated shapes for mass and angular projections as weighed by the fit proce- 
dure. As expected, the cos{9* ) distribution shows a signal shape which is 
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Fig. 4. Scatter plot of M^^ vs cos(0*^) for simulated events: background (top) and 
signal (bottom). 



much more peaked than the background to —1. The M.^^ distribution shows 
a well identified gaussian shape around the Ks mass for the signal, while the 
background has an asymmetric shape peaked at lower mass values. We count 
A^(27) = 600 ± 35 signal events out of 2280 events in the 2D plot. 



At the end of the analysis chain and as an independent check of the quality 
of the fit weights found, we show in Fig. 8. left the Xf/t distribution for data 
and MC as weighed by the scalar factors previously determined. A similar 
comparison is done also for the inclusive angular photon distribution (see 
Fig. 8. right). The latter distribution clearly indicates the need for a fiat angular 
dependence as expected by the uniform decay of a spin particle in two 
photons. 
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Fig. 5. Scatter plot of M^^ vs cos{0*^) for data after acceptance selection, QCAL 
veto and Xf/t ^^t- 

4 Branching ratio evaluation and systematics 



The branching ratio is evaluated with respect to the JMi{Ks —>■ 27r°) by count- 
ing the -f^i-crash tagged events with A^^ = 4 in the same sample as follows: 



BR{Ks 27) 



iV(27rO) eTOT{2l\tag) ^ 



where the total efficiencies have been evaluated by MC after Ks tagging. 
We have assumed that the ratio of trigger, event classification and tagging 
efficiencies between the two decays is one. The signal total efficiency is the 
product of the efficiencies for the acceptance selection, the QCAL cut and the 
cut. Each single one has been evaluated as conditioned efficiency. 

The acceptance selection efficiency for the signal, after tagging, is 



£seli2'j) 



^3.2±0.2,t,t±0.1 



^syst) /O, 



(2) 



this large efficiency is due to the angular coverage of the calorimeter, the low 
energy threshold used and the fiat angular distribution of the decay products. 
The systematic error assigned to this efficiency has been found by varying 
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Fig. 6. Distribution of the invariant mass of two photons, M^^ (MeV), for the 
Kl J J control sample with a neutral vertex in the fiducial region below the DC: 
data (top), MC samples (bottom). 

the data-MC correction curves of the cluster reconstruction efficiency. The 
efficiency of the QCAL cut, after tagging and acceptance, is found by MC to 
be ~ 99.96 % . However, on data we have to apply the corrections due to 
accidental losses described in sec. 3.3 to obtain: 

= X Cq = (96.45 ± O.OW ± 0.26syst.)%. (3) 



The MC efficiency of the applied Xfit cut is ~ (63.3 ±0.7)%. To evaluate 
a systematic error due to the data-MC difference in the Xfit scale, we have 
looked at the Kl 77 control sample by requiring loose Xf/t(KL) and 
angular cuts. A few percent contamination exists. By building Xfit cumulative 
distributions for data and MC and calculating their ratio at the applied cut 
value, a systematic error, A5^2/e^2 of -0.4% is assigned to this effect. 
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Fig. 7. Distributions of data (black points) after tagging, acceptance selection and 
Xfit '■ cos{6*^) (top), M^^ (bottom). The superimposed colored distributions 
are signal and background MC shapes. Solid line is the sum of the MC shapes after 
fitting. 

For the normalization we have counted Ks 27r° tagged events with N-y = 4. 
An efficiency of 

^sel 

(47) = (65.0 ± 0.2 

stat 

±0.1 

syst )% (4) 



is found by Monte Carlo. The systematics has been evaluated, as done for the 
signal, by varying the data-MC correction curves of the cluster reconstruction 
efficiency. After correcting for egg/ (47), a total number of (159.8 ±0.5) x 10^ 
Ks 27r° tagged events is obtained. Another systematic uncertainty related 
to accidental overlap of dafbkg clusters, shower fragmentation and merging 
of nearby clusters has been evaluated by repeating the measurement in an 
inclusive way and counting tagged events with 3, 4 and 5 photons. We get an 
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Fig. 8. Distributions of data (black points) after tagging, acceptance selection and 
Xfit • Xfit (left) and inclusive cos{9^) for the two photons in the event (right). 
The superimposed colored distributions are signal and background MC shapes. Solid 
line is the sum of MC shapes after fitting. 



Source 


+ABR/BR (%) 


-ABR/BR (%) 


Signal acceptance 


0.12 


0.12 


QCAL losses 


0.02 


0.26 


scale 




0.41 


Normalization sample 


0.15 


0.15 


QCAL TW change 


0.88 


0.44 


change 


0.44 


0.44 


2D-Fit binning 


0.88 


0.44 


MC Energy scale 




1.32 


Total 


1.33 


1.61 



Table 1 

Contributions to the total systematic error on the BR. The first three contributions 
have been evaluated directly as systematics related to the signal efficiency for a 
given cut. The fourth contribution regards the systematics on efficiency for the 
normalization sample. The last four contributions have been evaluated by repeating 
the BR measurement while varying analysis conditions or cuts. 



efficiency corrected counting of 159.5 x 10 events which agrees at ~ 2.5 per 
mil level with the number obtained with the exclusive counting. 
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Fig. 9. Variation of BR result as a function of the applied Xfit (left). Compar- 
ison of our measurement of JiIi{Ks — > 77) with the other existing measurements 
and 0{P^) predictions (right). 

For the BR{Ks 77) we use the latest PDG [10] value of BR{Ks 27r°) 
which is (30.69 ± 0.05)%. The systematics connected to the counting has been 
evaluated by repeating the analysis in different ways. We first tested the sta- 
bility of the branching ratio when modifying the width of the time window 
used for the QCAL veto or the applied value of the Xfit cut. In Fig. 9. left the 
BR chang function of the applied Xfit '^^^ shown. We then repeated 

the fit by applying the residual energy scale shift of +0.4% to the MC distri- 
butions and varied the bin size used in the 2D plot for fitting. In all cases, the 
maximum variation of the BR is used as systematic error and shown in Tab. 1. 
The sum in quadrature of all entries in the table is used as total systematic 
error. We obtain: 

BR{Ks 77) = (2.27 ± 0.13{stat.)toMisyst.)) x 10"^ (5) 



5 Conclusion 



With a sample of 1.62 fb~^ of e~^e~ collisions at ^/s ~ collected with KLOE 
at DA<I>NE, we have measured the BR{Ks — > 77) with a 5.6% statistical 
uncertainty and a ~ 1.5 % systematic error. We obtain a BR result which 
deviates by 2.9 cr's from the previous best precise determination, as shown in 
Fig. 9. right. Our measurement is also consistent, within errors, with O(P^) 
xPT predictions. 
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